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ABSTRACT 

The first deep blank-field 450/im map (lcr w 1.3 mJy) from the SCUBA-2 Cosmology 
Legacy Survey (S2CLS), conducted with the James Clerk Maxwell Telescope (JCMT) is 
presented. Our map covers 140arcmin 2 of the Cosmological Evolution Survey (COSMOS) 
field, in the footprint of the Hubble Space Telescope (HST) Cosmic Assembly Near-Infrared 
Deep Extragalactic Legacy Survey (CANDELS). Using 60 submillimetre galaxies (SMGs) 
detected at ^3.75ct, we evaluate the number counts of 450/im-selected galaxies with flux 
densities S450 > 5 mJy. The 8" JCMT beam and high sensitivity of SCUBA-2 now make it 
possible to directly resolve a larger fraction of the cosmic infrared background (CIB, peaking 
at A ~ 200/im) into the individual galaxies responsible for its emission than has previously 
been possible at this wavelength. At S450 > 5mJy we resolve (7.4 ± 0.7) x 10~ 2 MJy sr _1 
of the CIB at 450/im (equivalent to 16 ± 7% of the absolute brightness measured by the 
Cosmic Background Explorer at this wavelength) into point sources. A further ^40% of the 
CIB can be recovered through a statistical stack of 24/tm emitters in this field, indicating that 
the majority («60%) of the CIB at 450/im is emitted by galaxies with S450 > 2 mJy. The 
average redshift of 450/im emitters identified with an optical/near-infrared counterpart is es- 
timated to be (z) = 1.3, implying that the galaxies in the sample are in the ultraluminous 
class (Lir k 1.1 x 10 12 L Q ). If the galaxies contributing to the statistical stack lie at similar 
redshifts, then the majority of the CIB at 450/im is emitted by galaxies in the LIRG class with 
Lir > 3.6 x 10 n L o . 

Key words: galaxies: high-redshift, active, evolution, cosmology: observations, submillime- 
tre: galaxies 



1 INTRODUCTION 

Fifteen years have passed since the first 'submillimetre galaxies' 
(SMGs) were discovered (Smail, Ivison & Blain 1997; Barger et al. 
1998; Hughes et al. 1998), a high-redshift population (2 ~ 2 - 3, 
Chapman et al. 2005; Aretxaga et al. 2007; Wardlow et al. 2011) 
with ultraluminous (10 12 Lq) levels of bolometric emission, the 
bulk of which is emitted in the far-infrared (FIR) and redshifted to 
submillimetre wavelengths at z > 1. The power of submillimetre 
surveys for exploring the formation phase of massive galaxies was 
recognised before their discovery (e.g. Blain & Longair 1993; Dun- 
lop et al. 1994), and since their discovery, their importance as a cos- 
mologically significant population has been established by many 
studies (e.g. Smail et al. 2002; Dunlop et al. 2004; Ivison et al. 
2000, 2005, 2010; Coppin et al. 2008; Michalowski et al. 2010; 
Hainline et al. 201 1; Hickox et al. 2012; and see Dunlop et al. 201 1 
for a review). As such, SMGs provide challenging tests for models 
of galaxy formation, both in detailed 'zoomed' simulations as well 
as in cosmological theatres (Baugh et al. 2005; Dave et al. 2010). 
However, our view of the SMG population remains incomplete. 

In ground-based work, the majority of SMGs have - so-far - 
mainly been selected in the 850/iin or 1 mm atmospheric windows 
(e.g. Coppin et al. 2006; Weiss et al. 2009; Austermann et al. 2010; 
Scott et al. 2010), but this is far removed from the peak of the Cos- 
mic infrared background (CIB), which is at A ~ 200/im (Fixsen 
et al. 1998). The next available window closer to the CIB peak is 
at 450/im, but the transmission of this window is just at best 50% 
of the 850/iin window, making 450/im SMG surveys challenging 
from ground-based sites. Submillimetre surveys working closer to 
the CIB peak are essential if we are to identify the galaxies respon- 
sible for its emission; the 6*450 /Ssso colours of sources identified 
in the very deepest (lensing assisted) submillimetre surveys (Blain 
et al. 1999; Knudsen et al. 2008) suggests that these sources con- 
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tribute less than half of the CIB at 450/im (and therefore even less 
at the actual peak). 

The Balloon-borne Large Aperture submillimetre Telescope 
(BLAST, Pascale et al. 2008) made progress by conducting a low 
resolution submillimetre survey from the stratosphere at 250, 350 
and 500/rni (Pascale et al. 2008; Devlin et al. 2009; Glenn et 2010). 
This work was taken forward by the Herschel Space Observatory, 
which carries an instrument that images in the same wavelength 
ranges as BLAST (the Spectral and Photometric Imaging Receiver; 
SPIRE), and has mapped hundreds of square degrees of the sky at 
250-500/iin in a combination of panoramic and deep cosmological 
surveys (Eales et al. 2010, Oliver et al. 2010, 2012). However, the 
low resolution and high confusion limits of Herschel (FWHM~0.5' 
at 500/iin, a con w 7mJybeam~ 1 , Nguyen et al. 2010) limit the 
fraction of the CIB that can be directly resolved, with 15% resolved 
into individual galaxies at 250/iin and 6% at 500/iin (Clements et 
al. 2010; Glenn et al. 2010; Bethermin et al. 2012a). Thus, there 
remains work to be done in identifying the galaxies that emit the 
CIB, and thus finally complete the census of dust-obscured activity 
in the Universe and its role in galaxy evolution. 

Advances in submillimetre imaging technology are just now 
allowing us to take up the search once more, taking advantage of 
higher resolution possible with large terrestrial telescopes, and im- 
proved sensitivity and mapping capability in submillimetre detec- 
tor arrays. The SCUBA-2 camera is the state-of-the-art in submil- 
limetre wide-field instrumentation (Holland et al. 2006). The cam- 
era, now mounted on the 15 m James Clerk Maxwell Telescope 
(JCMT), consists of 5000 pixels in both 450/rni and 850/iin detec- 
tor arrays with an 8' x 8' field-of-view (16x that of its predeces- 
sor, SCUBA). The increase in pixel number is the reward of devel- 
opments in submillimetre detector technology; SCUBA-2 utilizes 
superconducting transition edge sensors (TES) to detect submil- 
limetre photons, with multiplexed superconducting quantum inter- 
ference device (SQUID) amplifiers handling read-out, analogous 
to an optical CCD. SCUBA-2 offers the capability to efficiently 
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Figure 1. (left) SCUBA-2 450^im signal-to-noise ratio map of the COSMOS/CANDELS field. The map has been scaled to emphasize the visibility of 60 
sources detected at >3.75cr significance (circled). The grey contours show the variation in the noise level, and are at 0-450 = 2, 3, 4, 5 mjy beam -1 (the solid 
angle bounded by the (7450 = 5mjybeam -1 contour is f2 fa 140arcmin 2 ). (right) Herschel-SVWE 500itm image of the same region, from the HerMES 
survey (Oliver et al. 2012). This map has been slightly smoothed with a Gaussian kernel to improve presentation. We show the limiting 0-459 = 5 mJy beam -1 
contour used for 450itm detection in the SCUBA-2 map and the position of the same galaxies in the left panel. This illustrates the ability of SCUBA-2 to 
push below the Herschel confusion limit at similar wavelengths, resolving confused emission into individual galaxies. 



map large (degree-scale) areas, and has the sensitivity to simulta- 
neously achieve deep (confusion limited) maps at both 450/im and 
850/xm. At 450/im, the resolution attainable with the JCMT is a 
factor ~5 x finer than the 500/iin resolution of Herschel, and the 
confusion limit is ~7 x fainter. 

Here we present results from early science observations of 
one of the seven components of the JCMT Legacy SurvejQ: the 
SCUBA-2 Cosmology Legacy Survey (S2CLSl| The goal of the 
S2CLS is to fully exploit SCUBA-2's mapping capabilities for the 
purpose of exploring the high redshift Universe. The S2CLS will 
cover several well-studied extragalactic 'legacy' fields, including 
the United Kingdom Infrared Deep Sky Survey Ultra Deep Sur- 
vey field (UDS), the Cosmological Evolution field (COSMOS), the 
Extended Groth Strip, and the Great Observatories Origins Deep 
Survey (North) fields. We present the deepest blank-field map at 
450/im yet produced (in the COSMOS field), and measure the flux 
distribution and abundance of the extragalactic sources revealed 
within it. In §2 we describe the observations and data reduction 
technique, in §3 we calculate the 450/im number counts and eval- 
uate the contribution to the CIB at 450/im. We briefly discuss and 
summarize our findings in §4 & §5. 



2 OBSERVATIONS AND DATA REDUCTION 

Observations were conducted in Band 1 weather conditions 
(T225GHz < 0.05) over 22 nights between 23 rd January and 20 th 
May 2012 totalling 50 hours of on-sky integration. The mapping 
centre of the SCUBA-2 COSMOS/CANDELS field is a = 10 h 



10 m 29.8 s , 5 = 02° 15' 01.6", chosen to be in the footprint of the 
Hubble Space Telescope CANDELS (Grogin et al. 2011; Koeke- 
moer et al. 20110 A standard 3 arcmin diameter 'daisy' mapping 
pattern was used, which keeps the pointing centre on one of the 
four SCUBA-2 sub-arrays at all times during exposure. 

2.1 Map making 

Individual 30min scans are reduced using the dynamic iterative 
map-maker of the SMURF package (Jenness et al. 2011; Chapin et 
al. 2012 in prep). Raw data are first flat-fielded using ramps brack- 
eting every science observation, scaling the data to units of pW. 
The signal recorded by each bolometer is then assumed to be a 
linear combination of: (a) a common mode signal dominated by at- 
mospheric water and ambient thermal emission; (b) the astronom- 
ical signal (attenuated by atmospheric extinction); and finally (c) 
a noise term, taken to be the combination of any additional sig- 
nal not accounted for by (a) and (b). The dynamic iterative map 
maker attempts to solve for these model components, refining the 
model until convergence is met, an acceptable tolerance has been 
reached, or a fixed number of iterations has been exhausted (in this 
case, 20). This culminates in time-streams for each bolometer that 
should contain only the astronomical signal, corrected for extinc- 
tion, plus noise. The signal from each bolometer's time stream is 
then re-gridded onto a map, according to the scan pattern, with the 
contribution to a given pixel weighted according to its time-domain 
variance (which is also used to estimate the \ 2 tolerance in the fit 
derived by the map maker). 



1 http ://w w w.j ach. hawaii . edu/JCMT/survey s 

2 http://www.jach.hawaii.edu/JCMT/surveys/Cosmology.html 
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Figure 2. Histogram of values in the SCUBA-2 450/tm signal-to-noise ra- 
tio map (Fig. 1), indicating the characteristic positive tail due to the pres- 
ence of real astronomical sources. The solid line is a Gaussian centred at 
zero with a width of a = 1, and the darker shaded histogram shows the 
histogram of pixel values in a map constructed by inverting a random 50% 
of the input scans; we use this for simulations of completeness, described 
in §2.3. Our detection limit is chosen to be a = 3.75, which yields a rea- 
sonably complete and reliable catalogue (see §2.3). Note that the 'real' map 
noise distribution is slightly wider than expected for pure Gaussian noise; 
this is due to slight 'ringing' around bright sources after convolution with 
the beam. 



The sky opacity at JCMT has been obtained by fitting ex- 
tinction models to hundreds of standard calibrators observed since 
the commissioning of SCUBA-2 (Dempsey et al. 2012). The op- 
tical depth in the 450/xm band was found to scale with the Cal- 
tech Submillimetre Observatory (CSO) 225 GHz optical depth as: 
T450 = 26.0(r 2 25-0.0196). Note that this scaling is slightly differ- 
ent from the original SCUBA relations (see Archibald et al. 2002; 
Dempsey et al. 2012). 

Filtering of the time-series is performed in the frequency do- 
main, with band-pass filters equivalent to angular scales of 2" < 
8 < 120" (i.e. frequencies of / = v/0, where v is the scan speed). 
The reduction also includes the usual filtering steps of spike re- 
moval (>10<7 deviations in a moving boxcar) and DC step correc- 
tions. Throughout the iterative map making process, bad bolome- 
ters (those significantly deviating from the model) are flagged and 
do not contribute to the final map. Maps from independent scans 
are co-added in an optimal stack using the variance of the data 
contributing to each pixel to weight spatially aligned pixels. Fi- 
nally, since we are interested in (generally faint) extragalactic point 
sources, we apply a beam matched filter to improve point source 
detectability, resulting in a map that is convolved with an estimate 
of the 450^im beam. The average exposure time over the nominal 
3 arcminute daisy mapping region (in practice there is usable data 
beyond this) is approximately lOksec per 2" x 2" pixel. 

2.2 Flux calibration 

The flux calibration of SCUBA-2 data has been examined by 
analysing all flux calibration observations since Summer 2011 un- 
til the date of observation. The derived beam-matched flux conver- 
sion factor (FCF) has been found to be reasonably stable over this 
period, and the average FCFs agree (within error) to those derived 



from the subset of standard calibrators observed on the nights of the 
observations presented here. Therefore we have adopted the canon- 
ical calibration of FCF450 = 540 ± 65 Jy beam -1 pW -1 here. A 
correction of ~ 10% is included in order to compensate for flux lost 
due to filtering in the blank-field map. This is estimated by inserting 
a bright Gaussian point source into the time stream of each obser- 
vation to measure the response of the model source to filtering. 

2.3 Maps and source detection 

We present the 450/im signal-to-noise ratio map of the COS- 
MOS/CANDELS field in Fig. 1. For comparison, we also show a 
Herschel SPIRE 500^im map of the same region to illustrate the 
gain in resolution that JCMT/SCUBA-2 offers at similar wave- 
lengthfl The 450/xm map has a radially varying sensitivity, which 
is nearly uniform over the central 3' (the nominal mapping area) 
and smoothly increases in the radial direction as the effective ex- 
posure time decreases for pixels at the edge of the scan pattern, 
which have fewer bolometers contributing to the accumulated ex- 
posure. The total area of the map considered for source extraction 
is 140arcmin 2 , where the rms noise is below 5mJybeam _1 . A 
histogram of pixel values in the 17450 ^ 5mJybeam _1 region is 
shown in Fig. 2. 

To identify extragalactic point sources, we search for pix- 
els in the (beam convolved) signal-to-noise ratio map with values 
>£ t hresh- If a P ea k is found, we record the peak-pixel sky co- 
ordinate, flux density and noise, mask-out a circular region equiv- 
alent to ~1.5x the size of the 8" beam at 450/im, reduce Ethresh 
by a small amount and then repeat the search. The floor value, be- 
low which we no longer trust the reality of 'detections' is chosen to 
be the signal-to-noise level at which the contamination rate due to 
false detections (expected from pure Gaussian noise) exceeds 5%, 
corresponding to a significance of <r~3.75. We detect 60 discrete 
point sources in this way, and these are identified in Fig. 1. Note 
that the map is far from confused, with an average source density 
equivalent to 6 x 10~ 3 beam -1 . We project that the confusion limit 
is at ~ 1 mjy beam - 1 

Completeness is estimated by injecting a noise model with ar- 
tificial point sources. To create maps with no astronomical sources 
but approximately the same noise properties of the real map, we 
generate jackknife realisations of the map where, in each fake map, 
a random half of individual scans have their signal inverted before 
co-addition (e.g. Weiss et al. 2009). Fig. 2 shows the equivalent his- 
togram of signal-to-noise ratio values in the jackknife map, which 
demonstrates the clean removal of astronomical sources, and the 
similarity with pure Gaussian noise. The recovery rate of sources 
as a function of input flux and local noise gives the completeness 
function: 10 5 fake sources in batches of 10 are inserted into the 
jackknife map, where each source selected from a uniform flux dis- 
tribution 1 < (5450/mJy) < 40. The 2-dimensional completeness 
function is shown in Fig. 3. 

In addition to the completeness correction, this technique al- 
lows us to estimate the noise-dependent flux boosting that occurs 
for sources with true fluxes close to the noise limit of the map, and 

4 The Herschel map was made from the Level 2.5 processed data products 
downloaded from the public Herschel Science Archive. The data were co- 
added with sky coverage used as an estimator for image noise level, and re- 
binned into the SCUBA-2 image reference frame, using nearest-neighbour 
sampling. The Icr noise level of this SPIRE map (including confusion) is 
6.2 mJy beam -1 



© 0000 RAS, MNRAS 000, 000-000 



S2CLS: 450 [im number counts and the CIB 5 



so we can construct an equivalent 'surface' in the noise-(measured) 
flux plane that can be use to de-boost the fluxes measured for point 
sources in the real map (Table 1). The typical de-boosting correc- 
tion is B <10%. Finally, the source detection algorithm is applied 
to each of the jackknife maps with no fake sources injected in or- 
der to evaluate the false positive rate, which we find to be 5%, in 
agreement with the false detection rate expected for a map of this 
size assuming fluctuations from pure Gaussian noise. 

A test for any bias in the recovery and correction of the source 
counts was performed in the following way. We populated the jack- 
knife maps with a model source count model (Bethermin et al. 
2012b) down to a flux limit of 5450 = 0.01 mjy. Sources were 
then extracted in exactly the same manner as the real data and com- 
pleteness and flux boosting corrected as described above and then 
compared to the input distribution. This procedure was repeated 
100 times and the average recovered source counts compared to 
the input model. The recovered differential and cumulative number 
counts were found to be consistent with the input number count 
realisations, indicating that our source detection and completeness 
corrections are not significantly biased. 



3 ANALYSIS 

3.1 Number counts of 450/im emitters 



In Table 1 and Fig. 4 we present the number counts at 450/im, cor- 
rected for flux boosting and incompleteness. The differential counts 
are well-described by a Schechter function: 




with S' = 10 mJy (fixed at a well-measured part of the flux distri- 
bution), N' = (490 ± 104) deg~ 2 and a = 3.0 ± 0.7. 

At flux densities above 20 mJy, the number counts from this 
survey are complemented by the equivalent measurements from 
Herschel surveys, which survey wider areas at 500/im to shallower 
depths, and so find the rarer, bright sources (nearby galaxies, ex- 
tremely luminous distant sources and gravitationally lensed galax- 
ies) that are not present in our map (Clements et al. 2010; Negrello 
et al. 2010). We focus on two Herschel surveys; HerMES, which 
has obtained confusion limited maps reaching a detection limit of 
S500 » 20mJy (Oliver et al. 2012) and the Herschel-PfiLAS sur- 
vey, which has mapped several hundreds of square degrees at a 
shallower depth (Eales et al. 2010). As Fig. 4 shows, our 450/im 
counts are in excellent agreement at ~20 mJy where the Herschel 
and SCUBA-2 CLS survey flux distributions meet. Below Her- 
schel's confusion limit the 500/im galaxy number counts have been 
inferred statistically, by both stacking (Bethermin et al. 2012b) and 
pixel fluctuation analyses (Glenn et al. 2010), again indicating con- 
sistency with the directly measured 450/mi number counts in ap- 
proximately the same flux regime. 

Recently, Chen et al. (2012) presented SCUBA-2 450/tm ob- 
servations of SMGs in the field of the lensing cluster A 370. The 
benefit of observing a lensing cluster is - provided a lens model is 
known - the ability to probe further down the luminosity function 
than would otherwise be possible for the same flux limit, with faint 
background sources boosted by the cluster potential. We compare 
the 'delensed' counts of 450/im emitters derived from 12 galaxies 
in the field of A 370 in Fig. 4, indicating broad agreement with our 
blank-field counts within the errors in the same flux range. After 
delensing, Chen et al. (2012) are able to probe slightly fainter than 




Figure 3. Completeness of the 450/tm catalogue as a function of local noise 
and input flux based on input-and-recovery simulations using jackknife re- 
alisations of the map noise. Modelling the completeness as a 2-dimensional 
function is required due to the radially varying sensitivity in the map (Fig. 
1). The same simulations allow us to estimate the difference between true 
(input) flux and recovered (i.e. observed) flux densities across the same pa- 
rameter space, and we use this information to correct the number counts 
accordingly. 

our catalogue, and the number counts at the 4.5 mJy level are also 
consistent with an extrapolation of our best fit Schechter function 
to the same limit. 

3.2 Resolving the 450/im background light 

What fraction of the CIB at 450/im have we resolved into galax- 
ies? The integrated flux density of point sources detected at 450/im 
(corrected for completeness) is J„(450/im) = (7.4 ± 0.7) x 
10 -2 MJy sr _1 . The absolute intensity of the CIB at 450/im mea- 
sured by C0B£-FIRAS is /„(450/xm) = 0.47 ± 0.19 MJy sr" 1 , 
thus we have directly resolved 16 ± 7% of the CIB at 450/im 
(the uncertainty is dominated by the COBii-FIRAS measurement; 
Fixsen et al. 1998). For comparison, the deepest Herschel surveys 
have directly resolved 5-6% of the CIB at 500/tm (Oliver et al. 
2010; Bethermin et al. 2012a). We show the integrated brightness 
of the 450/im emitters, relative to the absolute intensity of the CIB 
in Fig. 4. 

To measure the contribution to the CIB at 450/tm by galaxies 
not formally detected in the SCUBA-2 map, but which are known 
to be infrared-bright galaxies, we stack the map at the position of 
1600 galaxies selected from a catalogue generated from the Spitzer- 
COSMOS MIPS 24/im image of the same region (Sanders et al. 
2007). First, we remove point sources from the 450/xm map, using 
a point spread function (PSF) constructed by averaging the two di- 
mensional profiles of sources detected at >7<r. This PSF was then 
normalised to the flux of each individual source in our catalogue, 
and subtracted from the map. This yields a residual map where the 
only flux (in addition to that of noise) is contributed by sources not 
in our catalogue. The 450/im is then stacked at the position of the 
24/tm sources, averaging the flux with a weight equivalent to the 
inverse of the variance of the map at each position. 

The average 450/tm flux density of 24/im sources with mean 
24/jm flux (S24) = 0.19 mJy is (S450) = 2.2 ± 0.4mJy. 
The resulting contribution to the 450/im background is 0.20 ± 
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Figure 4. (left) Differential number counts of galaxies detected at 450/tm (error bars are derived from Poisson statistics). We compare the 450/tm counts to 
those measured recently by Chen et al. (2012) and by Herschel at 500/im (Bethermin et al. 2012a; Glenn et al. 2010; Oliver et al. 2010; Clements et al. 2010), 
and to the predictions of a numerical model of galaxy formation (Baugh et al. 2005; Lacey et al. 2008) and a parametric model of the evolution of the infrared 
luminosity density of the Universe (Bethermin et al. 2012b). We fit the measured 450/im counts with a Schechter function; the reason it fails to reproduce the 
shape of the 500/tm number counts (and those predicted by the models) at bright fluxes is because in this flux regime the counts have significant contributions 
from (a) bright local star-forming galaxies and (b) distant galaxies boosted to high observed flux by gravitational lensing; the SCUBA-2 map is too small 
to adequately sample these populations. Note - no 450/tm/500/im colour correction has been made to the 500/im data, (right) Integrated surface brightness 
of 450/im emitters relative to the CIB measured by COBE-FIRAS at 450/jm (Fixsen et al. 1998). The directly measured number counts are well-fitted by a 
Schechter function, the extrapolation of which agrees with the CIB derived from a stack of 24/im-emitting galaxies not individually detected in the SCUBA-2 
map. Thus, we directly resolve 16±7% of the CIB measured by COB£-FIRAS into galaxies, with an additional !%40% contributed by 24/im-emitting galaxies 
not formally detected at 450/im. We project that 100% of the CIB at 450/im is recovered at a flux density of 5450 > 0.4 mjy. 



0.04 MJy sr -1 (the uncertainty is a/y/~N, with a the standard de- 
viation in the stack and N the sample size). A simple simulation 
was performed to test whether the stacking methodology described 
above produces unbiased estimates of the submillimetre flux. The 
residual flux map was inverted (by multiplying by —1) and sim- 
ulated sources were inserted using the derived PSF as a model, 
with the input fluxes of the fake sources set to S450 = 10524 up 
to a maximum of S450 = 5mJy. The positions were set to the 
24/im catalogue positions, rotated 90 degrees about the map cen- 
tre, thus preserving clustering information. The stacking procedure 
was then repeated as for the real catalogue. The mean input flux 
was S450 = 1.8 mjy per source, and the recovered stacked flux was 
S450 = 1.0 ± 0.5 mJy. The recovered flux is slightly low compared 
to the input flux at the 1.5cr level, however this does not affect our 
conclusions, given the uncertainties in the 450/im flux calibration 
and the absolute measured value of the CIB at 450/im. 

Excluding those detected as bright point sources, the 24/im- 
selected galaxies contribute (2.0 ± 0.4) x 10 _1 MJy sr -1 , or 42 ± 
19% of the CIB at 450/im. Therefore, in addition to the directly 
detected sources, in total we can account for 58 ± 20% of the 
CIB at 450/im using the SCUBA-2 map. Note that our stacked 
value is in good agreement with the background derived from a 
stack of 24/tm-emitters in Herschel 500/tm images (Bethermin et 
al. 2012b), and is also consistent with the intensity expected from 
an extrapolation of a Schechter function fit to the directly measured 
number counts (Fig. 4). 



Table 1. Number counts of 450/tm-selected galaxies. N indicates the raw 
number of galaxies in each bin (8S45Q = 5 mJy), and the completeness (C) 
and de-boosting (B) corrections represent the mean corrections for galaxies 
in each bin (note that each galaxy is de-boosted individually, with the cor- 
rection increasing for lower flux densities). Uncertainties in the differential 
counts are the lcr confidence range assuming Poisson statistics (Gehrels 
1986). 



5 45 o N dN/dS N(>S') a (C) h (B) c 

(mjy) (mJy" 1 deg~ 2 ) (deg~ 2 ) 



8.0 


41 


343.0±*j;» 


231 3 4+ 339 - 7 


1.6 


1.1 


13.0 


13 


QQ K + 32.3 


508 5+ 172 '4 


1.1 


1.1 


18.0 


1 


20.8±^ 8 




1.0 


1.0 


23.0 


2 




52.otS:S 


1.0 


1.0 



a 5' corresponds to the lower edge of the bin, i.e. (S450 — 2.5) mjy 
b average completeness correction applied 
c average flux de-boosting correction applied 



4 DISCUSSION 

We compare our results to the phenomenological model of 
Bethermin et al. (2012b), who use a 'backwards evolution' param- 
eterisation of the the infrared luminosity density (as traced by dusty 
star-forming galaxies; see also Lagache et al. 2004). The Bethermin 
et al. (2012b) model assumes that the star formation modes of 
galaxies can be either described as 'main sequence' (i.e. SFR scales 
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with stellar mass) or 'starburst', with spectral energy distributions 
defined by the latest stellar synthesis template libraries. The evo- 
lution of the luminosity functions of these two populations inte- 
grated over cosmic history provides good fits to the observed num- 
ber counts of galaxies at 24, 70, 100, 160, 250, 350, 500, 850, 
1 100/im and 1 .4 GHz (as well as integrated observables such as the 
evolution of the volume averaged star formation rate and cosmic 
infrared background). Here we confirm that the number counts of 
450/im emitters predicted by the model is also in good agreement 
with the measured 450/im number counts in the flux range probed 
by our SCUBA-2 survey. 

We also compare the measured counts to the GALFORM semi- 
analytic model of galaxy formation (Cole et al. 2000; Baugh 2006; 
Lacey et al. 2008; Almeida et al. 2011). This prescription predicts 
the formation and evolution of galaxies within the A CDM model 
of structure formation (Springel et al. 2005), and includes the key 
physics of the galaxy formation (and evolution) process: radiative 
cooling of gas within the dark matter halos, quiescent (by which 
we mean non-burst driven) star formation in the resultant discs, 
mergers, chemical enrichment of the stellar populations and inter- 
galactic medium and feedback from supernovae and active galactic 
nuclei. As Fig. 4 shows, the numerical model slightly over-predicts 
the abundance of 450/im emitters in the flux range probed. Never- 
theless, the reasonable agreement between the shape of the counts 
predicted by GALFORM and the data is encouraging for models of 
galaxy formation that aim to reproduce the full range of emission 
processes of galaxies at long wavelengths. 

The 8 arcsec resolution of the 450/im SCUBA-2 map allows 
us to accurately identify the optical/near-infrared counterparts of 
the SMGs, and we have identified the most likely counterpart to 
the majority of 450/im sources in our sample (I. G. Roseboom et al. 
2012 in prep). The wealth of legacy data available in the COSMOS 
field then provides the means to estimate the redshift distribution 
of the population. We have used 13 bands of optical/near-infrared 
photometry, including CFHT ugri, Subaru SuprimeCam z', VISTA 
YJHK, HSTF125W and F160W and Spitzer IRAC [3.6] and [4.5] 
to evaluate the photometric redshifts of all the identified galaxies 
(the typical la uncertainty based on the confidence level of the 
template fit is 8z = 0.16). The redshift distribution is shown in 
Fig. 5, indicating that the majority of our sample lie at z < 3, with 
a mean redshift of (2) = 1.3 (a full analysis of the source iden- 
tification and redshift distribution is to be presented in Roseboom 
et al. 2012 in prep). This is a clear indication that the 450/im se- 
lection is probing a lower redshift population than previous 850/im 
selected samples, which have typical redshifts of (z) — 2.2 (e.g. 
Chapman et al. 2005; Wardlow et al. 2010). The shape of the red- 
shift distributions predicted both by the phenomenological model 
and numerical model described above (for galaxies at the same flux 
limit) are also in good agreement with the measured distribution; 
both models predict little contribution from galaxies at z > 3 (al- 
though a high redshift 'tail' is present in both models). 

Assuming the directly detected sources representing 16±7% 
of the CIB at 450/im are star-forming galaxies at (z) = 1.3, then 
their total (rest-frame 8-1000/im) luminosities are in the ultralu- 
minous class, L m « 1-1 X 10 12 L Q (Chary & Elbaz 2001). If the 
galaxies contributing to the 24/iin stack described in §3.3 lie in 
the same redshift range as the directly detected galaxies, then the 
majority (60%) of the CIB at 450/im is emitted by galaxies with 
Lm. > 3.6 x 10 11 Lq. This is broadly consistent with the picture 
that at z ~ 1, the star formation rate budget of the Universe is dom- 
inated by galaxies in the LIRG class, with star formation rates of 
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Figure 5. Redshift distribution of 450/tm-selected SMGs in our sample (not 
corrected for completeness), derived from 13-band photometric redshift es- 
timates (§3.4). The average redshift of the sample is (2)= 1 .3, and the vast 
majority of the 450/tm-selected SMGs lie at z < 3. For comparison, the 
average redshift of SMGs selected at 850/mi is z 2.2 (Wardlow et al. 
2010), indicating the efficacy at which the 450/tm selection samples a pop- 
ulation of SMGs at lower redshift and therefore an important complement to 
any census of the dusty Universe. We compare the shape of the redshift dis- 
tribution to the models of Bethermin et al. (2012b) and Lacey et al. (2008) 
shown in Fig. 4 for galaxies with S450 > 5 mJy (we have area-normalised 
both model distributions, since the observed redshift distribution contains 
no completeness correction). The average redshift and shape of both model 
distributions is in good agreement with observations, suggesting that - at 
this flux limit - there is little contribution from galaxies at z > 3. 

order 10 M Q yr" 1 (Dole et al. 2006; Rodighiero et al. 2010; Mag- 
nellietal. 2011). 

An extrapolation of the Schechter function fit to the directly 
measured number counts (which agrees well with the background 
at 5450 ~ 2mJy, derived from the stack of 24/im sources), im- 
plies that 100% of the CIB at 450/im should be recovered at 0.1 < 
S450 < 1.4 mJy (the range accounting for the la uncertainty of 
the absolute measured background; Fixsen et al. 1998), close to the 
SCUBA-2 confusion limit. If the galaxies responsible for this emis- 
sion are at similar redshifts to the current 450/im sample (but below 
the sensitivity of the map and not contributing to the 24/im stack), 
then the majority of the remaining ~40% of the CIB at 450/im is 
likely to be emitted by galaxies with Lm < 1-3 x 1O 11 L0, im- 
plying galaxies star formation rates of a few tens of Solar masses 
per year. However, we cannot as yet rule out what fraction of the 
remaining CIB light might be emitted by faint 450/im emitters at 
higher redshifts; note that a galaxy with S450 ~ 2mJy at z > 2 
has a typical luminosity of Lir > 5.5 x 10 11 Lq (Chary & Elbaz 
2001), again indicating the importance of LIRG-class galaxies in 
the cosmic infrared budget. Characterizing the high redshift tail of 
the 450/im population is an important next step. 



5 SUMMARY 

The SCUBA-2 camera on the 15 m JCMT repesents the state-of- 
the-art in panoramic submillimetre imaging, and has recently begun 
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scientific observations in earnest. In this paper we have presented 
results from the first deep, blank-field cosmological map at 450/im 
(o"45o = 1.3 mJy); part of the SCUBA-2 Cosmology Legacy Sur- 
vey, the largest of the seven JCMT Legacy Surveys. Using a 450 /im 
map of the well-studied extragalactic COSMOS/CANDELS field, 
we have 

(i) made the first unbiased, blank-field determination of the 
number counts of galaxies at 450/im, at a flux density limit of 
S450 > 5mJy. This probes below the confusion limit of Her- 
schel, complementing the number counts measured at fluxes above 
20 mJy over wider areas in major Herschel submillimetre surveys; 

(ii) measured the contribution of these galaxies to the cosmic in- 
frared background at 450/im: we resolve 16% of the CIB into indi- 
vidual galaxies. The ability of SCUBA-2 to 'pin-point' the galaxies 
responsible for the emission of the CIB is a critical step in under- 
standing the properties of the galaxies that are forming the majority 
of stars in the Universe at this epoch; 

(iii) an additional ~40% of the CIB can be recovered in the 
SCUBA-2 map by stacking Spitzer MlPS-detected 24/^m emitters. 
Using this analysis we estimate that the majority («60%) of the 
CIB at 450/im is emitted by galaxies with S450 > 2 mJy; 

(iv) a preliminary analysis of the redshift distribution of the 
450/im emitters (based on high-quality photometric redshifts avail- 
able for this field) imply that the typical redshift of galaxies with 
S450 > 5mJy is (z) = 1.3, with the majority lying at 2 < 3. 
The typical luminosity of galaxies in our sample are estimated to 
be in the ultraluminous class, with Lir > 1O 12 L0. If the galax- 
ies contributing to the statistical stack of 24/im emitters described 
above are at a similar redshift, then we project that the majority 
of the CIB at 450/im is emitted by 'LIRG' class galaxies with 
Lir > 1.3 x 10 n L Q . 

These are the first results of the S2CLS. The final goal of the survey 
will be to map a quarter of a square degree to (7450 = 1.2 mJy, and 
a wider, ten square degree area to ogso = 1.5 mJy, yielding >10 4 
SMGs with which to (a) determine the submillimetre luminosity 
function and its evolution over cosmic time; (b) search for the 
rarest, most luminous SMGs at high-z; (c) resolve the 450/im back- 
ground; (d) accurately measure the clustering properties of SMGs 
and determine their typical environments (including rare protoclus- 
ters) and (e) build-up the large samples required to properly re- 
late SMGs to other star-forming (ultraviolet/near-infrared selected) 
populations at z ~ 2 and thus gain further insight into SMGs' role 
in the overall history of galaxy evolution. 
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